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An Antisymmetrized Molecular Dynamics model is used to explore the collision of 114 Cd projec- 
tiles with 92 Mo target nuclei at E/A=50 MeV over a broad range in impact parameter. The atomic 
number (Z), velocity, and emission pattern of the reaction products are examined as a function 
of the impact parameter and the cluster recognition time. The non-central collisions are found to 
be essentially binary in character resulting in the formation of an excited projectile-like fragment 
(PLF*) and target-like fragment (TLF*). The decay of these fragments occurs on a short timescale, 
100<t<300 fm/c. The average excitation energy deduced for the PLF* and TLF* 'saturates for 
mid-central collisions, 3.5<b<6 fm, with its magnitude depending on the cluster recognition time. 
For short cluster recognition times (t=150 fm/c), an average excitation energy as high as ~6 MeV 
is predicted. Short timescale emission leads to a loss of initial correlations and results in features 
such as an anisotropic emission pattern of both IMFs and alpha particles emitted from the PLF* 
and TLF* in peripheral collisions. 

PACS numbers: PACS number(s): 25.70.Mn 
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I. INTRODUCTION 



Collision of two heavy-ions at intermediate energies 
can result in the production of a multi-particle final state 
PI 0, Q • These multi-particle final states have been ex- 
perimentally characterized by a wide variety of signals 
including fragment multiplicity EL EL size distributions 
E M , emission timescales @, Q, H, H, scaling behavior 
|lCt If ll [f^ | and the attained excitation energy 0] • F° r 
large fragment multiplicity, within a thermodynamic ap- 
proach, such multi-fragment states have been interpreted 
as a transition of the finite nuclear system from a liquid to 
a gaseous phase US- EU 113 ■ Recent work has investi- 
gated the robustness of this conclusion by examining the 
influence of the surface, through the density dependence 
of the entropy, on the stability of the nuclear droplet 
against fragmentation 0,^1- All these approaches how- 
ever focus on the thermodynamic stability of the system. 
In reality, the decaying system is formed by the colli- 
sion dynamics which may not equilibrate all degrees of 
freedom equally [2fJ, ■ ^ n order to understand both 
the formation and decay of excited nuclear systems in- 
volved in the collision process, microscopic app roaches 
have also been followed [H |2l H |H Hi . In or- 
der to make direct comparison with experimental data 
such microscopic models typically utilize a a two-stage 
approach. In the first phase, a dynamical model is used 
to describe the collision dynamics. Clusters produced 
in this phase are subsequently de-excited by a statistical 
model. Such a two-stage approach typically views the 
statistical decay stage as decoupled from the dynamical 
stage that preceded it. In the present work we exam- 
ine the validity of such a de-coupled hybrid approach. 
Specifically, we utilize a microscopic model, the Anti- 



symmetrized Molecular Dynamics model, to investigate 
how the collision proceeds on short timescales and how 
the reaction characteristics evolve with impact parame- 
ter. In addition, we examine whether initial correlations, 
existing at short times, survive the decay stage and how 
they are manifested in final distributions. 



II. DESCRIPTION OF THE AMD MODEL 

To describe the dynamical stage of intermediate en- 
ergy heavy-ion collisions, we utilize the antisymmetrized 
molecular dynamics (AMD) model jH 1^3, . For 
the present work, we use the same version of AMD as 
Ref. [3(j which has been used to describe the multifrag- 
mentation reaction of the central Xe + Sn collisions at 50 
MeV/nucleon. 

The description of the dynamics of fragmentation is, 
in principle, a very complicated quantum many-body 
problem. In the exact solution of the many-body time- 
dependent Schrodinger equation, the intermediate and 
final states should be very complicated states contain- 
ing a huge number of reaction channels corresponding to 
different fragmentation configurations. The AMD model 
respects the existence of channels, while it neglects some 
of the interference among them. Namely, the total many- 
body wave function |\&(t)) is approximated by 
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w(Z,t)dZ, (f) 



where each channel wave function |$(Z)) is parametrized 
by a set of parameters Z , and w(Z, t) is the time- 
dependent probability of each channel. 
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In AMD, we choose the Slater determinant of Gaussian 
wave packets as the channel wave function 



exp{-i>(r; - Zj/V^ 2 }xa 3 'W 
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where x Qi are the spin-isospin states with oti — p j 
,p J.,n t, or n j. Thus, the many-body state 
is parametrized by a set of complex variables Z = 
{Zi}i=i,... t A, where A is the number of nucleons in the 
system. The width parameter, v = 0.16 fm -2 , is treated 
as a constant parameter common to all the wave packets. 
If we ignore the antisymmetrization effect, the real part 
of Zj corresponds to the position centroid and the imag- 
inary part corresponds to the momentum centroid. This 
choice of channel wave functions is suitable for fragmen- 
tation reactions, where each single particle wave function 
should be localized within a fragment. 

Instead of directly considering the probability w(Z,t) 
in Eq. QJ, we solve a stochastic equation of motion for 
the wave packet centroids which may be symbolically 
written as 



d 
lit 



Z l = {Z t , H}pb + (NN coll) + AZ 4 (t) +/x (Z* , «')• (3) 



The first term {Zj,7i}pB is the deterministic term de- 
rived from the time-dependent variational principle with 
an assumed effective interaction. The Gogny interaction 
|32| is used in the present work. The second term repre- 
sents the effect of the stochastic two-nucleon collision pro- 
cess, where a parametrization of the energy-dependent 
in-medium cross section is adopted. The two-nucleon 
collision cross-section used is the same as in Ref. |3fT | 
namely, 



a(E,p) = minf cr LM (E,p), 



100 mb 



1 + £/(200 MeV) 



(4) 



The collisions are performed with the "physical nucleon 
coordinates" that take account of the antisymmetriza- 
tion effects, and then the Pauli blocking in the final state 
is automatically introduced [28l |2^ |. The third term 
AZj(t) is a stochastic fluctuation term that has been 
introduced in order to respect the change of the width 
and shape of the single particle distribution [23, [33|> 13 • 
In other words, the combination {Zj, H}pb + AZj(i) ap- 
proximately reproduces the prediction by mean field the- 
ories (for a short time period) for the ensemble-averaged 
single-particle distribution, while each nucleon is local- 
ized in phase space for each channel. The term AZ,*(t) 
is calculated practically by solving the Vlasov equation 
(for a short time period) with the same effective interac- 
tion as for the term {Zj, H}pb- In the present version of 
AMD |3(j], the property of the fluctuation AZj(t) is cho- 
sen in such a way that the coherent single particle motion 
in the mean field is respected for some time interval un- 
til the nucleon collides another nucleon. The last term 
fi (Zj, TC) is a dissipation term related to the fluctuation 
term AZj(t). The dissipation term is necessary in order 



to restore the conservation of energy that is violated by 
the fluctuation term. The coefficient /i is given by the 
condition of energy conservation. However, the form of 
this term is somehow arbitrary. We shift the variables 
Z to the direction of the gradient of the energy expec- 
tation value TL under the constraints of conserved quan- 
tities (the center-of-mass variables and the total angular 
momentum) and global one-body quantities (monopole 
and quadrupole moments in coordinate and momentum 
spaces). A complete formulation of AMD can be found 
inRefs. mill. 



The statistical decay of relatively small primary frag- 
ments (Z < 20) is calculated by using the code [35| based 
on the sequential binary decay model by Piihlhofer |36| . 
The code employed in the present work also takes ac- 
count of the emission of composite particles not only in 
their ground states but also in their excited states with 
the excitation energy E* < 40 MeV. The experimental 
information is incorporated for known levels of A < 28 
nuclei, while the Fermi-gas level density is assumed other- 
wise. For the statistical decay of large primary fragments 
(Z > 20), the decay code GEMINI [33 is employed. In 
considering the decay of the fragments, both the exci- 
tation energy and decay probabilities are calculated for 
spherical fragments independent of the true shape of the 
fragments induced by the reaction dynamics. The effect 
of n-p asymmetry, excitation energy, and deformation on 
the nuclear level density are not considered in the de- 
cay. Introduction of a deformation dependence of the 
nuclear level density, and in particular the treatment of 
the continuum, results in a significant modification of the 
emission rate for fragments that are weakly bound or at 
high excitation |38| ]. 

The system we have chosen to study is 114 Cd + 92 Mo 
at E/A = 50 MeV which can be considered represen- 
tative of symmetric heavy-ion collisions in this energy 
domain. We sampled all impact parameters, b, in the 
interval 0<b<b max with a triangular distribution. The 
maximum impact parameter b max bad a value of 12 
fm. The touching sphere configuration distance, given 
by R=1.2*(A P 1 / 3 +A T 1 / 3 ), is equal to 11.2 fm. The pro- 
jectile and target were therefore placed at an initial dis- 
tance of 13 fm for b>6.5 fm and 9.8 fm for b<6.5 fm. 
For a given collision, the fate of the colliding system was 
followed until 300 fm/c. At regular intervals, the posi- 
tions and momenta of all nucleons in the system were 
recorded. At a selected time (typically 300 fm/c), which 
we designate the cluster recognition time, the nucleon 
distributions are subjected to a cluster recognition al- 
gorithm based on the distance between nucleons. The 
nucleons and clusters that result from cluster recognition 
are subsequently propagated along Coulomb trajectories 
and allowed to statistically decay. The identity and mo- 
menta of the final reaction products are recorded for sub- 
sequent analysis. In order to examine the predictions of 
this model in a statistically significant manner, we have 
amassed «25,000 collisions. The calculations were per- 
formed on a 646 CPU parallel processor system of which 
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each CPU was either a PowerPC or Power3+. A single 
collision for this reaction required 12 to 24 CPU-hours 
on a node depending on the impact parameter. 



III. GENERAL REACTION 
CHARACTERISTICS 

Depicted in Fig. is the density distribution of nucle- 
ons in R-space as a function of time for a mid-peripheral 
(b=7.79 fm) and mid-central (b=5.15) collision. The ini- 
tial moment in time (t=0) is taken as the near touching 
configuration of the projectile-target system previously 
described, with the projectile approaching the target nu- 
cleus from the negative z direction. As the di-nuclear 
system rotates, the initial dumb-bell shape of the two 
touching nuclei shown in the top panel evolves. While in 
contact, the two nuclei exchange mass, charge, and en- 
ergy, governed by nucleon-nucleon scattering within the 
mean field. 

For the presented event with b=7.79 fm, one observes 
that two large nuclei emerge from the collision at t=90 
fm/c revealing the intrinsically binary nature of the col- 
lision. In this case, at longer times the elongation of 
the target-like fragment (left) leads to its breakup into 
multiple intermediate size nuclei. In the case of the mid- 
central collision with b=5.15 fm however, the situation 
is more difficult to discern. At t=90 fm/c, it is un- 
clear whether the system is disassembling into two or 
three large pieces. What is apparent is that as the two 
nuclei separate from each other, one observes that the 
density distributions reflect the nuclear interaction be- 
tween the projectile and target nuclei through the forma- 
tion of transiently deformed nuclei. These non-spherical 
geometries persist up to 300 fm/c for different cluster 
sizes. Moreover, for both events presented clusters seem 
to emerge on a relatively fast timescale, t«90 fm/c. This 
early production of clusters indicates that the timescale 
of the shape/density fluctuations responsible for cluster 
formation operate on this timescale. It should be noted 
that a considerable fraction, though not all, of this early 
stage cluster emission is located between the two large 
fragments that emerge from the collision. The evolution 
of the density distributions presented in Fig. ^ can also 
be viewed from the context of semi-classical colliding liq- 
uid drops. Formation of the transiently extended nuclear 
system by the collision dynamics involves the generation 
of a considerable amount of "surface" nuclear material as 
compared to "bulk" nuclear material. In comparison to 
the original system comprised of the projectile and tar- 
get nuclei, the multi-fragment final state with multiple 
clusters requires the formation of a significant amount 
of surface - an energetically unfavorable change. Thus, 
once the surface-to-volume ratio has been increased by 
the collision dynamics, the energy cost of the system re- 
organizing to the multi-fragment final state is consider- 
ably reduced. 

We examine the characteristics of the system imme- 
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FIG. 1: Contour diagram depicting the nucleon density dis- 
tribution in spatial coordinates as a function of time for the 
reaction 114 Cd + 92 Mo at E/A = 50 MeV. The positive z 
direction corresponds to the direction of the projectile. The 
columns correspond to different impact parameters, b=7.79 
fm (left) and b=5.15 fm (right). 



diately following this dynamical stage of the collision. 
The products of the reaction at this stage are designated 
the "primary" products which statistically de-excite to 
form the final reaction products which we also refer to 
as the "secondary" products. For a large ensemble of 
events we examine the evolution of both primary and 



4 



secondary distributions with impact parameter, velocity 
dissipation, and cluster recognition time. 




V„ (cm/ns) 



(cm/ns) 



FIG. 2: (Color online) Two dimensional diagram of the cor- 
relation between the atomic number and parallel velocity of 
particles at t=300 fm/c for different impact parameters. The 
arrows correspond to the projectile and target velocities. The 
color scale indicates the yield on a logarithmic scale. 

An overview of the collisions studied is presented in 
Fig. |21 where the correlation between the atomic number 
and parallel velocity (in the center-of-mass frame) of par- 
ticles at t=300 fm/c is examined. For the most peripheral 
collisions (8.5<b<10 fm) two peaks located at Z«47 and 
«39 are clearly evident. These peaks correspond to the 
excited projectile-like (PLF*) and target-like (TLF*) nu- 
clei respectively and are relatively narrow distributions 
in velocity centered at Vy =3.8 and -4.7 cm/ns. Also ev- 
ident is copious production of neutrons (Z=0), hydro- 
gen, and helium nuclei. Smaller in yield, are clusters 
with Z>3 and atomic number less than that of the PLF* 
and TLF*. This pattern, dominated by the survival of 
the PLF* and TLF* for a peripheral collision, reflects a 
primarily binary nature. For mid-peripheral and mid- 
central collisions, a similar pattern is observed indicating 
that in this impact parameter range as well a PLF* and 
TLF* survive the dynamical phase, hence these impact 
parameters are also essentially binary in character. With 
increasing centrality (Vplf*) decreases and (Vtlf*) in- 
creases reflecting an increase in the velocity damping. At 
the same time, the width of the PLF* and TLF* velocity 
distributions increases indicating the growth of fluctu- 
ations. In addition, with increasing centrality the av- 
erage atomic number of the PLF* and TLF* decreases 
while the yield of clusters with 3<Z<15 increases. For 
simplicity, we designate the highest Z cluster with a ve- 




FIG. 3: Dependence of the primary (solid) and secondary 
(dashed) Z distributions on impact parameter. The differen- 
tial yield dN/dZ has been normalized by the total number of 
events for each impact parameter interval. 



locity larger (smaller) than the center-of-mass velocity 
as the PLF* (TLF*). For b<4 fm the decrease in the 
average Z of the PLF* combined with the width of the 
distribution, lead to an operational definition of interme- 
diate mass fragment, namely IMF: 3<Z<10. Particles 
with Z<10, manifest broad velocity distributions for the 
most central collisions. Examination of the most periph- 
eral collisions reveals a clear pattern of how the velocity 
distribution evolves with the atomic number (Z) of the 
fragment. Neutrons and hydrogen nuclei in particular 
have velocity distributions that are centered on veloci- 
ties between those of the PLF* and TLF*. In contrast, 
for nuclei with 3<Z<15 the velocity distribution while 
broad, clearly has a bimodal nature with each of the two 
peaks centered close to the PLF* and TLF* velocities. 
This bimodal character is also observed for helium nu- 
clei although the distributions are broader. These overall 
patterns manifested for the most peripheral collisions are 
also observed for more central collisions. 

Depicted in Fig. [3] is the dependence of the primary 
and secondary Z distributions on impact parameter. The 
Z distribution of particles at t=300 fm/c is the primary 
distribution and is represented as the solid histogram. 
Following Coulomb propagation and statistical decay of 
the excited primary reaction products, the Z distribution 
of secondary particles is represented by the dashed his- 
togram. The latter distribution includes both primary 
fragments that did not decay, as well as the decay prod- 
ucts of excited primary fragments. All distributions have 
been normalized to the total number of events for each 
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impact parameter range and therefore represent the av- 
erage multiplicities. As may be expected from the trends 
in Fig. |3 the charge distribution for the most peripheral 
collisions, b>10 fm, is largely dominated by two peaks at 
Z=42 and 48 which correspond to the TLF* and PLF*. 
In Fig. |2| we therefore focus on smaller impact parame- 
ters, b<10 fm. In panel a) the Z distributions integrated 
over impact parameter up to 10 fm are presented. As ex- 
pected, the yield for neutrons, hydrogen, and helium is 
large in the primary distribution (solid histogram). A 
large yield is also observed for 3<Z<10. Evident for 
Z>30 is a slight double peak in the primary distribu- 
tion attributable to the presence of the PLF* and TLF* . 
This double peak structure is eliminated by secondary 
decay as it is not evident in the dashed histogram. To 
separate the PLF* from the TLF* , as well as to crudely 
separate their decay products, we further select particles 
with the condition V|| >0. The resulting primary distri- 
bution shown in panel b) manifests only a single peak at 
large Z, which is located at Z=47. As observed in panel 
a) the yield of the Z distribution for 3<Z<30 is similar 
for both the primary and secondary particles. 

We examine the dependence of the Z distribution on 
impact parameter for V||>0 in Fig. [3]>f. For 8.5<b<10 
fm, Fig. |3J:), the primary Z distribution is 'V-shaped', 
reminiscent of the 'U-shape' observed for asymmetric fis- 
sion. The minimum yield observed near Z«20 is deep 
in comparison to the yield at lower and higher Z in- 
dicating that asymmetric splits are strongly preferred 
over symmetric splits. It is striking that the multiplic- 
ity for Z=3-6 is approximately the same as that of Z«47 
(the PLF*). The yield ratio for Z=3-6 over Z=45-47 is 
0.31/0.37~0.84, indicating a process or processes result- 
ing in copious production of light IMFs. This similarity 
in the yield of the light IMF and the PLF* can, for exam- 
ple, be understood as the asymmetric binary decay of a 
precursor PLF*. Such a perspective is supported by ex- 
perimental observation. For peripheral collisions of two 
heavy-ions at intermediate energies, the phenomenon of 
dynamical fission is well characterized [33, tHi] • This 
dynamical fission has been associated with the deforma- 
tion of the PLF* induced by the collision process. The 
defining characteristics of this process are the aligned 
asymmetric binary decay of the PLF* and large rela- 
tive velocities between the two produced fragments. On 
general grounds one expects that this dynamical process 
should depend sensitively on both the induced deforma- 
tion and the excitation of the PLF* 42]. It is impor- 
tant to observe that the shape of this primary distribu- 
tion largely survives the process of secondary decay. The 
main difference between the primary and secondary dis- 
tributions is that the high Z peak is shifted to lower Z 
and increases in width. For 7<b<8.5 fm, Fig. |3Ji), the 
shape of the primary distribution is better described as a 
'U-shape'. In contrast to the previous impact parameter 
interval, the minimum located at Z«20 is shallow. This 
decrease in the depth of the minimum can be associated 
with the increase in the probability of symmetric binary 



splits relative to asymmetric binary splits. This change 
of the Z distribution with decreasing impact parameter 
can be related to an increase in the excitation energy of 
the PLF*. In this impact parameter interval, the yield 
for Z=3-6 is significantly larger than that for Z«42. The 
ratio of the yield of Z=3-6 over the yield of Z=41-43 is 
0.66/0. 19~3. 47, a change by a factor of «4 as the impact 
parameter decreases from 8.5<b<10 fm to 7<b<8.5 fm. 
This increase in the ratio is due to both an increase in 
the IMF yield by a factor of 2 and a decrease in the yield 
in the vicinity of the PLF* peak. The latter decrease 
reflects the increasing width of the peak in the Z distri- 
bution attributable to the PLF* with decreasing impact 
parameter. Following secondary decay the 'U-shape' is 
somewhat less pronounced. For yet more central col- 
lisions, a 'U-shape' distribution is not observed even for 
the primary distribution. In panel e) no clear bump is ob- 
served at large Z, indicating the decreased likelihood that 
a high Z PLF* survives to the cluster recognition time of 
t=300 fm/c. For the most central collisions shown, b< 4 
fm, the primary Z distribution is exponential over a large 
range in Z. This exponential behavior of the yield is sup- 
pressed for Z>30 due to the finite size (atomic number) 
of the system. The secondary Z distribution for central 
collisions also exhibits an exponential character for Z>3, 
although the onset of the finite size effects is observed 
at Z=20. The main effect of secondary decay on the Z 
distribution, for all impact parameters is to significantly 
enhance the yield of neutrons, hydrogen, and helium nu- 
clei, while decreasing the maximum Z observed, namely 
the atomic number of the PLF* and TLF*. 

We have investigated whether the similarity of the 
yield for Z=3-6 and Zp^p* for peripheral collisions is 
an indication that the two largest fragments forward 
of the center-of-mass originate from a common parent. 
Displayed in the two-dimensional diagrams of Fig. 0] is 
the joint probability of observing the largest and sec- 
ond largest fragments both with Vy>0. For reference the 
dashed lines correspond to Ztot = %plf* + T^i^Uggest 
= 45, 40, 35, 30 and 25. The distribution at t=300 fm/c 
is presented in the left column while the distribution fol- 
lowing secondary decay is shown in the right column. In 
the case of primary fragments (left column), for b>4 fm, 
a clear anti-correlation is observed between the atomic 
number of the largest and second largest fragment. In 
order to examine the average behavior of the two dimen- 
sional distribution, we also indicate as solid and open 
circles the (^•z nd Uggest) for a given Zp^p. and (Zplf») 
for a given 'Zi2 nd biggest- Strong correlation of (Zp£p«) and 
('Z'2 nd biggest} is evidenced by the near overlap of the open 
and closed circles over an extended range. Divergence of 
the symbols indicates that either the two fragments do 
not originate from a common parent or that finite size 
effects strongly influence the observed correlation. For 
8.5<b<10 fm, Ztot is almost constant over the range of 
Zp_LP* with a value of « 45. This value corresponds to an 
average loss of three charges from the incident Cd nucleus 
with the observed anti-correlation signaling a conserva- 
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FIG. 4: (Color online) Correlation between the atomic num- 
ber of the fragments with the largest and second largest 
atomic number that have velocities larger than the center-of- 
mass velocity. The solid circles represent the average charge of 
the second biggest fragment for a given charge of the biggest 
fragment. The open circles represent the average charge of 
the biggest fragment for a given charge of the second biggest 
fragment. 



tion of charge between the largest and second largest frag- 
ment. This anti-correlation signals that both fragments 
do on average originate from a common parent fragment. 
For a second largest fragment with Z=6, the (Zplf*) is 
«40 consistent with dynamical breakup 0, l43| . While 
mid-peripheral (mid-central) collisions exhibit an anti- 
correlation, Ztot changes from «45 (45) at high Zplf* 
to w40 (35) at low Zplf*- This change in Ztot might 
indicate that the PLF* splits into three or more pieces 
or simply reflect the changing size of the parent fragment 
over the finite impact parameter interval considered. For 
b<4 fm, the average atomic number of the largest and 
second largest fragment are closer, (Z)—20 and 9 respec- 
tively, as one might expect for a binary decay at high 
excitation, consistent with the increased probability of 
symmetric splits indicated by the Z distribution in Fig. [3] 
The general trends observed for the primary fragments 
are also evident following decay in the charge correlation 
of secondary fragments (right column). Similar charge 
correlation patterns have been experimentally observed 
|40| indicating a transition from asymmetric splits to- 
ward those in which all asymmetries are populated. The 
total charge of the two fragments after secondary decay 




b (fm) b (fm) 

FIG. 5: (Color online) Dependence of the (Z),(V), and E*/A 
of the PLF* and TLF* on impact parameter. Error bars 
indicate the standard deviation of the distribution. In the 
bottom panels the solid symbols denote the ((E* / A)xlf*) as 
a function of b. 



is typically reduced by 5 to 10 charges as compared to 
the total charge of the primary fragments. Following sec- 
ondary decay, the emission of Z=4-6 results in a marked 
horizontal line in the charge correlation. This feature in 
the charge correlation has also been experimentally ob- 
served and has been previously attributed to dynamical 
fission [ZjJ. 



IV. VELOCITY DISSIPATION OF THE PLF* 
AND ITS EXCITATION 

A more quantitative picture of the evolution of the 
general properties of the PLF* and TLF* with impact 
parameter is displayed in Fig. [S] For the most periph- 
eral collisions studied, b>10 fm, the (Zplf*) is «48, the 
atomic number of the projectile. The (Zplf*) decreases 
smoothly with decreasing impact parameter until b~3- 
4 fm. For smaller impact parameters, (Zplf*) shows no 
dependence on impact parameter and has a value of «19. 
For b<10 fm, (Zplf), namely the average atomic num- 
ber following decay, is approximately 4-9 units less than 
(Zplp*) and exhibits the same impact parameter depen- 
dence as (Zplf*)- K should be noted that the largest 
difference between (Zplf*) and (Zplf) is observed for 
mid-peripheral collisions with an impact parameter «8 
fm. The average center-of-mass velocity of the PLF*, 
(VpLF*), also exhibits a smooth dependence on impact 
parameter, decreasing monotonically from (Vplf*) ~ 4.3 
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cm/ns for the most peripheral collisions to «2.5 cm/ns 
for b=3 fm. For more central collisions (Vplf*) only 
shows a weak dependence on impact parameter. With 
increasing centrality the width of the velocity damping 
distribution (indicated by the error bars) increases sig- 
nificantly, indicating the growth of fluctuations. 

The predicted velocity damping of the PLF* evident 
in the middle panel is associated with a corresponding 
increase in the excitation of the PLF* as shown in the 
bottom panel of Fig. [SJ Such an association between ve- 
locity damping and excitation has been experimentally 
observed Q. While the average E*/A of the PLF* 
rapidly increases with impact parameter for peripheral 
collisions, it saturates at «4 MeV by b=6 fm. The trends 
observed for the PLF* are also observed for the TLF* as 
depicted in the right column of Fig. It is interesting 
to note that the (E*/A) for small impact parameters at- 
tained for both the PLF* and TLF* is the same despite 
the smaller size of the TLF* (Z«15) as compared to the 
PLF* (Z«19). This difference of «20-25% in Z corre- 
sponds to a similar difference in A (see Fig. EJ) . Equal 
partition of E* would thus result in a larger (E*/A) for 
the TLF* as compared to the PLF*. An (E*/A)=4 MeV 
for the PLF* would correspond to an (E*/A)=5 MeV for 
the TLF*. The similarity of (E*/A) for both the PLF* 
and TLF* is indicative that the degree to which thermal- 
ization is achieved is large. For the most peripheral colli- 
sions, b~12 fm, the non-zero value of the ((E*/A)tlf*) 
and ((E*/A)pii?*) is due in part to the mismatch be- 
tween the binding energy of the projectile and target in 
AMD and their real binding energies. This error typi- 
cally ranges from 0.2 to 0.3 MeV. Additional excitation 
may occur due to the mean field or Coulomb interaction. 

The composition of the excited PLF* and TLF* that 
subsequently undergoes decay is indicated in Fig. |SJ In 
this figure the (N/Z) of both the PLF* and TLF* are ex- 
amined as a function of impact parameter. For b>6 fm, 
the (N/Z) of both PLF* and TLF* evolves essentially 
linearly with impact parameter from the initial values 
of 1.375 and 1.19 for the projectile and target respec- 
tively. Over this range of impact parameter, this behav- 
ior could be interpreted as equilibration of N/Z. How- 
ever, the change in (N/Z) is larger for the PLF* as com- 
pared to the TLF* by a factor of two. This difference 
reflects the fact that exchange between the PLF* and 
TLF* is not the only process occurring thus complicat- 
ing the interpretation of the change in N/Z in terms of 
equilibration. For more central collisions, the (N/Z) re- 
mains essentially constant having saturated at a value of 
sal. 24-1. 26. The similarity of the average N/Z value for 
the PLF* and TLF* could be interpreted as equilibra- 
tion of this degree-of-freedom. If this is indeed the case, 
it is interesting to note that for b«4 fm, this equilibra- 
tion is already achieved. For comparison the (N/Z) of 
the system is «1.29. The slightly lower N/Z asymptotic 
value for central collisions as compared to the N/Z of the 
system suggests either a preferential emission of free neu- 
trons or the production of neutron-rich fragments in the 
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FIG. 6: Top panel: Dependence of the (N/Z) of the PLF* and 
TLF* on impact parameter. Dotted lines indicate the initial 
N/Z of projectile and target nuclei. Bottom panel: Fraction, 
f, of nucleons found in the PLF* (solid symbols) or TLF* 
(open symbols) that originate from the projectile. 



dynamical stage. 

We examine the degree to which mixing occurs in the 
lower panel of Fig. [S] In this figure the dependence of f, 
the fraction of nucleons in the PLF* or TLF* that orig- 
inate from the projectile, on impact parameter is pre- 
sented. It is interesting to note that for b>6 fm, the 
region in which (N/Z) changed linearly with b, the frac- 
tion of nucleons in the PLF* that were originally in the 
projectile is large, f>0.9. Only for smaller impact param- 
eters does the degree of mixing of projectile and target 
nucleons become larger. Thus, the large change in (N/Z) 
does not require large mixing of the projectile and target 
nucleons. It is instructive to note that the quantity f, ap- 
pears to saturate for b<2 fm with a maximum of «35% 
of the PLF* nucleons originating from the target. For 
the TLF*, in the case of small impact parameters, the 
degree of mixing is similar. It has been experimentally 
demonstrated that for mid-peripheral collisions the N/Z 
degree of freedom does not reach equilibrium |45| . How- 
ever, the present result indicate that N/Z equilibrium is 
attained for mid-central collisions, despite the incomplete 
mixing of the projectile and target nucleons. This result 
is of significance to future work with radioactive beams, 
indicating the degree to which the N/Z exotic projectile 
can be excited while only modestly pertubing its N/Z. 

To probe the origin of the saturation in E*/A of the 
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FIG. 7: Dependence of the (E*/A) of the PLF* on impact 
parameter for different cluster recognition times. 



PLF* and TLF* for b<6 fm observed in Fig. El we have 
investigated the influence of our choice of cluster recog- 
nition time on the excitation energy of the PLF*. We 
have chosen to recognize the clusters at t= 150, 180, 
210, 240, 270, and 300 fm/c and compare the depen- 
dence of excitation energy on impact parameter for the 
different cluster recognition times. As evident in Fig. 
while for peripheral collisions the average excitation en- 
ergy is fairly independent of the choice of cluster recog- 
nition time, with decreasing impact parameter the aver- 
age excitation energy deduced depends significantly on 
the choice of cluster recognition time. For different clus- 
ter recognition times one also observes that the onset 
of the saturation in excitation energy occurs at different 
impact parameter. For t=300 fm/c the onset of the sat- 
uration occurs at b«6 fm («25 % of the cross-section) 
while for t=150 fm/c, the onset occurs at b«4 fm («10 
% of the cross-section). The events associated with the 
highest excitation attainable therefore correspond to a 
significant fraction of the cross-section. For central col- 
lisions, the excitation attained is higher the earlier one 
recognizes the clusters. For an early cluster recognition 
time, t— 150 fm/c, a maximum value of (E*/A)«6 MeV 
is attained in comparison to (E*/A)ss4 MeV for t=300 
fm/c. Both the trend and magnitude of (E*/A) is consis- 
tent with the AMD calculations for a more asymmetric 
system j22|. This rapid decrease in (E*/A) is indicative 
of rapid cooling of the PLF* . As one may imagine, the 
choice of a cluster recognition time less than 150 fm/c 
becomes increasingly problematic due to both the con- 
ceptual, as well as practical, problem of distinguishing 
clusters during the high density phase of the collision. 

The dependence of (Vplf*) and (E*/A) on impact pa- 
rameter suggests a direct correlation between these two 
quantities. The correlation between these two quanti- 



ties as a function of both impact parameter and clus- 
ter recognition time is examined in Fig. |SJ For periph- 
eral collisions, 8.5<b<10 fm (leftmost column), at early 
cluster recognition times, e.g. t=150 fm/c (uppermost 
panel) , a narrow anti-correlated distribution is observed, 
namely there is a strong dependence of the PLF* 's excita- 
tion, (E* / A)plf* , on its velocity, Vplf*- To more eas- 
ily examine the correlation between the two quantities, 
the centroid in E*/A for each bin in Vplf* is indicated 
by the symbol. The significant slope of ((E* / A)plf*) 
with respect to Vplf* indicates the strong correlation 
between PLF* excitation and velocity damping. With 
increasing cluster recognition time, the strong correla- 
tion between Vplf* and (E* / A)plf* persists although 
the width of the distribution increases. 

For a fixed cluster recognition time, one observes that 
with decreasing impact parameter, the dependence of 
((E* /A)plf*} (symbols) on velocity becomes flatter in- 
dicating a weakening dependence on average. The two 
dimensional distributions also become broader with in- 
creasing centrality indicating the growth of fluctuations 
that attenuate the intrinsic correlation between excita- 
tion energy and PLF* velocity. Examination of the most 
central collisions studied (b<4 fm) shows that while a 
modest dependence between E*/A and Vplf* exists at 
t=150 fm/c, for longer cluster recognition times effec- 
tively no dependence of the PLF* excitation energy on 
its velocity is observed. At t=300 fm/c, ((E*/A) PLF .) 
does not exhibit any dependence on Vplf*- This atten- 
uation of the correlation between excitation energy and 
velocity of the PLF* with increasing cluster recognition 
time is also observed at intermediate impact parameters. 

To examine the influence of the cluster recognition 
time on the most peripheral collisions in a more quan- 
titative manner, we compare in the top panel of Fig. 
the dependence of the average excitation energy as a 
function of Vplf* for 8.5<b<10 fm for different clus- 
ter recognition times. For low to modest velocity damp- 
ing, i.e. Vpl_f*>3.25 cm/ns, the anti-correlation be- 
tween the average excitation energy and Vplf* is in- 
dependent of the cluster recognition time. For more 
damped collisions, however, one does observe a difference 
between the calculated average excitation energy for dif- 
ferent cluster recognition times. For V plf*—^ cm/ns, 
((E*/A) FiF .)=2.15 MeV while for V PL f*=3.5 cm/ns, 
((E*/A) PiF »)=3.3 MeV. This average excitation of 3.3 
MeV is associated with a velocity damping from beam 
velocity of 0.86 cm/ns. 

In the lower panel of Fig. the dependence of the av- 
erage excitation energy as a function of Vplf* for more 
central collisions 4<b<7 fm is shown. In contrast to the 
more peripheral collisions just discussed, for all values of 
velocity damping, the average excitation of the PLF* de- 
pends on the cluster recognition time. Even for the small- 
est velocity damping (Vplf*> 3.5 cm/ns) a minimum 
excitation energy of «3.5 MeV is observed for all cluster 
recognition times. Cluster recognition times less than 210 
fm/c manifest an essentially linear dependence of E*/A 
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FIG. 8: (Color online) Correlation between the E*/A of the PLF* and its velocity for different impact parameters and cluster 
recognition times. The symbols indicate the average E*/A as a function of Vplf* ■ The arrows correspond to the beam velocity. 



on V plf* while longer cluster recognition times (t>240 
fm/c) exhibit a significantly non-linear dependence. By 
t=300 fm/c, (E*/A) is essentially independent of Vplf* ■ 
The dependence of some of the average properties of 
the PLF* on both impact parameter and cluster recog- 
nition time are summarized in Fig. I1UI In the top panel, 
the average atomic number of the PLF*, (Zplf*), is dis- 
played as a function of cluster recognition time for dif- 
ferent impact parameters. For the most peripheral colli- 
sions, 8.5<b<10 fm, and the shortest cluster recognition 
times, (ZpLF*)&4:7, just below Zbeam=^& as indicated 
by the dashed line. Longer cluster recognition times re- 
sult in a slight decrease in (Zplf*) to a value of w 44 
at t=300 fm/c. This reduction in (Zplf*) corresponds 
to the emission of charge on a short timescale. For more 



central collisions a similar behavior is observed although 
the magnitude of the charge emitted on a short timescale 
is larger. 

In the middle panel of Fig. ^3 the trend of (Vplf*) 
with cluster recognition time and impact parameter is 
presented. For 8.5<b<10 fm, essentially no change is ob- 
served in (Vplf*) as the cluster recognition time changes 
from t=150 fm/c to 300 fm/c. For mid-central and cen- 
tral collisions, a small increase in (Vplf*) is discernible 
as the cluster recognition time increases. This slight in- 
crease is attributable to the Coulomb re-acceleration of 
the PLF* following the collision combined with recoil ef- 
fects due to predominantly backward emission of parti- 
cles on a short timescale. 

The dependence of ((E*/A)plf*) on cluster recogni- 
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FIG. 9: Average E*/A of the PLF* as a function of its velocity 
for different impact parameters and cluster recognition times. 
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on cluster recognition time for different impact parameters. 
The dashed line represents the projectile atomic number (ve- 
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tion time is depicted in the bottom panel of Fig. ITUI 
for different impact. As previously noted in Fig. for 
8.5<b<10 fm the cluster recognition time has only a 
weak influence on ((E* / A)plf»)- Longer cluster recog- 
nition times lead to slightly lower ((E* /A)plf*), 3.1 
MeV for t=150 fm/c as compared to 2.8 MeV at 300 
fm/c. More central collisions, however, manifest a more 
marked dependence. As apparent in Fig. [7| for b<4 fm 
((E* /A) plf*) reaches a value of 6 MeV for the shortest 
cluster recognition times, while at longer cluster recogni- 
tion times {(E*/A)plf*) is only w4 MeV. This decrease 
in excitation energy is rapid with most of the decrease 
occurring from t=150-240 fm/c. This rapid decrease in 
the excitation energy of the PLF* is directly related to 
the emission of particles over this time interval. As the 
excited PLF* rapidly emits charged particles between 
t=150-240 fm/c its atomic number decreases while its 
velocity remains relatively constant. Consequently, the 
correlation between ((E* /A)plf*) and Vplf* observed 
in Fig. 13 is poor for central collisions and long cluster 
recognition times. Thus, a proper description of this de- 
excitation of the highly excited PLF* requires modeling 
the statistical decay of the deformed PLF*. 

An interesting consequence of this rapid emission from 
the PLF* (and TLF*) is the amelioration of the correla- 
tion between the excitation of the PLF* and TLF*. Dis- 
played in Fig. 1111 is the two-dimensional distribution of 



PLF* and TLF* excitation energies for different impact 
parameters and cluster recognition times. For all impact 
parameters shown, the distribution is broad with the cen- 
troid for each (E* / 'A)p£f* bin indicated by the symbol. 
For 8.5<b<10 fm, a slight positive correlation between 
((E*/A)tlf*) and (E*/A) PLF * is evident. Examination 
of the correlation between the total excitation, E*, of 
the PLF* and TLF* reveals an independence indicating 
that the observed correlation between (E* / A)plf* and 
{{W / A)tlf*) is principally due to a correlation between 
Aplf* and Atlf* for the most peripheral collisions. In 
contrast, for the most central collisions, b<4 fm, an anti- 
correlation between ((E*/A) T lf*) and (E* /A)plf* is 
evident. Both E* and A of the PLF* and TLF* man- 
ifest the same anti-correlation. 



V. EMITTED PARTICLES 

As evident from Fig. [T] as the PLF* and TLF* sep- 
arate, clusters are produced. This fragment production 
as already demonstrated can occur on relatively short 
time scale impacting the Z, velocity, and (E*/A) of the 
PLF* and TLF*. In order to characterize this fast emis- 
sion process in more detail, we examine the multiplicity 
of fragments produced as a function of impact param- 
eter in Fig. Displayed in Fig. 112b ) is the average 
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FIG. 11: (Color online) Correlation between E*/A of the PLF* and TLF* for different impact parameters and cluster recognition 
times. The symbols correspond to {(E*/A)tlf* ) for a given (E*/A)p_lf*. 



multiplicity of fragments, Z>3, at t=300 fm/c (solid cir- 
cles). One observes that this multiplicity increases with 
decreasing impact parameter and saturates for b~3 fm. 
For the most peripheral collisions the average multiplic- 
ity is 2, corresponding to the existence of the only PLF* 
and TLF*. The average fragment multiplicity reaches a 
value of 3 at b~8-9 fm. For this impact parameter inter- 
val, on average, one fragment is produced in coincidence 
with the PLF* and TLF* . This result is consistent with 
the asymmetric split of the PLF* deduced from the Z dis- 
tribution (Figs. |21 and 0J). For the most central collisions, 
b<3 fm, the average fragments multiplicity is constant 
and is «8. Following secondary decay (solid triangles) 
the fragment multiplicity is reduced slightly due to the 
decay of fragments into particles with Z<2. For b>7 



fm the effect of secondary decay on the fragment multi- 
plicity is negligible while for the most central collisions 
the average multiplicity decreases from 8 to 6.5. The 
increased excitation energy associated with more central 
collisions is no doubt responsible for this increased impor- 
tance of secondary decay. The multiplicities predicted in 
the present calculation are compared to those from QMD 
calculations for the system Xe+Sn ^(|. Although both 
systems were simulated for the same incident energy of 
50 MeV/nucleon, the Xe+Sn system is «20% larger in A 
and ~15% larger in Z than the present system. The mul- 
tiplicity deduced by QMD (open crosses) is larger that 
the ones of the present work at all impact parameters. 
Given the difference in the system size, the difference be- 
tween the multiplicities for b<6 fm may be reasonable. 
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FIG. 12: (Color online) Multiplicity of fragments as a function 
of impact parameter. Panel a: Average multiplicity of Z> 
3 at t=300 fm/c (solid circles), after secondary decay (solid 
triangles) and for QMD (open crosses). The QMD results 
are extracted from |4q|. Panel b: Multiplicity of Z=3-10. 
Panel c: Average multiplicity of Z=3-10 for different cluster 
recognition times. 



The most notable feature of this comparison between the 
two models is the behavior for peripheral collisions, b>8 
fm. The fragment multiplicities predicted by AMD ap- 
pear to be more realistic than those predicted by QMD. 
This difference may be due to spurious decay of the pro- 
jectile and target in QMD due to the poor description of 
the ground state properties in that model. 

The multiplicity distribution of IMFs (3<Z<10) as a 
function of impact parameter is presented in Fig. 112b ) for 
t=300fm/c. While the distribution is narrow for the most 
peripheral collisions, its width rapidly increases with de- 
creasing impact parameter. The average IMF multiplic- 
ity, indicated by the solid circles, evolves from for the 
most peripheral to «6 for b=3 fm. At an impact parame- 
ter of «9 fm, the average IMF multiplicity reaches a value 
of sal, consistent with Fig. 112b). The average IMF mul- 
tiplicity is pretty insensitive to secondary decay as indi- 
cated by the triangles. Comparison between the fragment 
multiplicity, Fig. THk). and IMF multiplicity, Fig. H2b). 
indicates that even for the most central collisions two 
fragments with a Z>10 are present at t=300fm/c repre- 
senting a PLF* and TLF* with approximately 20-25 % 
of the original projectile and target atomic number. This 
result contradicts the physical picture of a single source 



often assumed for central collisions. 

Displayed in Fig. 112b ) is the average IMF multiplicity 
dependence on b for different cluster recognition times. 
The average IMF multiplicity increases with increasing 
cluster recognition time for all impact parameters . The 
largest increases are evident for the shortest times, t<240 
fm/c. For all impact parameters, the IMF multiplicity 
increases by a factor of 2 to 3 between t=150 fm/c and 
t=240 fm/c. After t=240 fm/c, the IMF production rate 
is reduced with an increase of 20-40 % of the IMF mul- 
tiplicity between t=240 fm/c and t=300 fm/c. 

The average multiplicity of light charged particles is 
examined in Fig. 1131 as a function of both impact param- 
eter and cluster recognition time. In the left hand column 
of Fig. EH one observes a monotonic increase of the neu- 
tron and proton average multiplicities with decreasing 
impact parameter both at t=150 fm/c and t= 300 fm/c. 
At t=150 fm/c a slight saturation in the both the neu- 
tron and proton multiplicities is observed for the most 
central collision with maximum average multiplicities of 
11.5 and 7 attained. A later cluster recognition time of 
t=300 fm/c results in approximately a 50% increase in 
the multiplicities with the saturation of the multiplici- 
ties for central collisions being slightly more evident. For 
this longer cluster recognition time, the average multi- 
plicities associated with central collisions are 19.5 and 12 
for neutrons and protons respectively. Following sequen- 
tial decay (t=oo), one observes a significant increase in 
the average multiplicities and a pronounced saturation in 
the case of the neutrons. This saturation suggests that 
the total neutron multiplicity, in particular, while provid- 
ing impact parameter selectivity for peripheral collisions 
is a poor selector of more central collisions. Moreover, 
attempting to select central collisions with the neutron 
multiplicity would on the basis of the cross-section be 
weighted towards mid-central collisions. This result ex- 
plains the experimental observation of the persistence of 
binary collisions associated with large neutron multiplic- 
ity |4jJ . Moreover, these mid-central collisions are associ- 
ated with the highest average excitation energy attained 
as presented in Fig. |SJ 

In the case of the light cluster (d,t, and a particle) 
multiplicities a couple of points are noteworthy. For 
short cluster recognition time (t=150 fm/c) the average 
multiplicity of deuterons is relatively linear over the en- 
tire impact parameter range and reaches a value of « 1.4 
for the most central collisions. Alpha particles, in the 
case of peripheral collisions manifest similar multiplici- 
ties, however the average multiplicity of a particles sat- 
urates for b<6 fm. Tritons exhibit lower multiplicities 
than both deuterons and alpha particles for all impact 
parameters. For longer cluster recognition time, t=300 
fm/c, the deuteron and triton multiplicities remain es- 
sentially unchanged as compared to t=150 fm/c. In con- 
trast, the a particle multiplicity increases significantly. 
It is interesting to note that the maximum a multiplicity 
is not associated with central collisions but rather with 
b«7 fm. From this we conclude that significant a produc- 
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FIG. 13: Multiplicity of light particles as a function of impact 
parameter for different cluster recognition times. The multi- 
plicities in the right column have been scaled by the factors 
indicated. 



tion/emission, but not deuteron or triton emission occurs 
on the timescale commensurate with the separation time 
ofthePLF* andTLF* (150 fm/c<t<300 fm/c). Follow- 
ing secondary decay (t=oo) all multiplicities increase sig- 
nificantly. Moreover, only for peripheral collisions, b>8 
fm, does the average multiplicity of light clusters depend 
significantly on impact parameter. 

We examine the emission pattern for IMFs (3<Z<10) 
in Fig.^Jboth at t=300 fm/c and at t=oo as a function 
of impact parameter. In examining the most peripheral 
collisions for t=300 fm/c, we observe two major compo- 
nents which are shifted with respect to the velocity of 
the PLF* and TLF* as represented by the arrows in the 
figure. In addition a minor component is visible centered 
at the velocity of the center-of-mass i.e. Vm =0. This 
emission pattern is consistent with anisotropic emission 
in the frame of the PLF* and TLF*. The most likely ori- 
gin if the observed backward enhancement, i.e. towards 
mid-rapidity, is the asymmetry of the collision process it- 
self. With increasing centrality, one observes an increase 
in this backward yield, as well as an increase in the yield 
of the mid- velocity component. For b<7 fm, this mid- 
velocity yield becomes considerable. For the most cen- 
tral collisions, the distinct bimodal character evident in 
more peripheral collisions is replaced by a broad distribu- 
tion. The impact of Coulomb propagation and secondary 
decay is shown in the right column of Fig. 1141 In con- 



trast to the broad distributions observed at t=300 fm/c, 
the emission pattern following Coulomb propagation to 
infinite PLF-TLF separation and secondary decay (right 
column) , reveals a pattern of two semi-circles centered on 
the PLF* and TLF* velocities. Such an emission pattern 
reflects both the Coulomb focusing in the field of the sep- 
arating PLF* and TLF* , as well as emission of IMFs from 
the de-exciting PLF* and TLF* . For the most peripheral 
collisions one observes two distinct Coulomb circles. It is 
important to note that the intensity pattern along each of 
these Coulomb circles is not constant but exhibits a sig- 
nificant backward enhancement indicating a memory of 
the initial angular asymmetry. With decreasing impact 
parameter, the center of these Coulomb circles shifts to- 
ward the center-of-mass and increasingly overlap as the 
velocity of the PLF* and TLF* decrease. The Coulomb 
circles also become less distinct with increasing central- 
ity reflecting both increased excitation of the system and 
nucleon-nucleon scattering. 

The parallel velocity distributions of the PLF*, TLF*, 
and IMFs and their decay products are shown in Fig. 1151 
as a function of impact parameter. The velocity distribu- 
tions of the PLF* and TLF* (left column) are presented 
for reference (shaded histogram) . For clarity these latter 
distributions have been scaled relative to the IMF distri- 
butions by the factors indicated. In the case of 8.5<b<10 
fm, the PLF* and TLF* manifest gaussian-like veloc- 
ity distributions that are relatively narrow and slightly 
damped from the beam velocity. With decreasing impact 
parameter, these two distributions move closer in veloc- 
ity, i.e. exhibit increased damping, and become broader. 
The parallel velocity distributions of the PLF and TLF 
(right column) follow the same general trends as those of 
the PLF* and TLF* . The widths of the secondary large 
fragments are typically 10-40 % larger than that of the 
PLF* and TLF*. 

For the most peripheral collisions, the IMF veloc- 
ity distribution (solid histogram) is bimodal with the 
most probable values of this two peaked distribution dis- 
placed toward the center-of-mass velocity as compared to 
the PLF* and TLF* velocities, clearly establishing the 
qualitative trend first observed in Fig. In addition 
to the two gaussian yields attributable to the emission 
from the PLF* and TLF* an additional IMF component, 
smaller in magnitude, is observed. As previously noted in 
Fig. O this additional component has an average veloc- 
ity roughly centered at the center-of-mass velocity. For 
7<b<8.5 fm, the relative magnitude of the mid-velocity 
contribution is increased. With increasing centrality, the 
shape of the IMF velocity distribution evolves toward a 
flat distribution reflecting increased fragment production 
at mid-velocity. 

We have fit the predicted parallel velocity distributions 
shown with two gaussians representing the emission from 
the PLF* and TLF* . The result is depicted as the dashed 
histogram in Fig. The fit parameters for the PLF* 
and TLF* emission are presented in Table [|] With in- 
creasing centrality the centroid of the PLF* velocity dis- 
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FIG. 14: (Color online) Invariant cross-section for IMFs (3< Z<10) in the COM frame. The arrows indicate the average parallel 
velocity of the PLF* and TLF*. The color scale indicates the yield on a logarithmic scale. The vertical scale in the bottom 
two panels has been scaled by a factor of two as compared to the other panels. 



TABLE I: PLF* and TLF* average parallel velocity, and fit 
parameters for the two gaussian fit of the IMF Vii distribu- 
tions at t=300 fm/c. The deduced quantities, (V\\) and a are 
expressed in cm/ns. 





PLF* 


IMF(PLF*) 


TLF* 


IMF(TLF') 


b (fm) 


<V||> 


<V||> 


a 


(V||> 


<V||> a 


8.5<b<10 


3.76 


2.59 


1.19 


-4.65 


-4.01 0.84 


7<b<8.5 


3.51 


2.93 


1.01 


-4.31 


-4.03 0.86 


4<b<7 


2.99 


3.02 


0.87 


-3.57 


-3.86 0.84 


b<4 


2.26 


2.55 


0.94 


-2.48 


-2.81 1.25 



TABLE II: Average multiplicity of 4<Z<Z PLF - , Ztlf* for 
the PLF* , TLF* and mid- velocity (MV) components at t=300 
fm/c. The relative fraction of the mid-velocity component to 
the total emission of the PLF* and TLF* is also indicated. 



b (fm) 


PLF* 


TLF* 


MV 


P(MV) 


3.5<b<10 


0.09 


0.09 


0.06 


0.26 


7<b<8.5 


0.63 


0.63 


0.81 


0.39 


4<b<7 


0.99 


0.94 


1.53 


0.44 


b<4 


1.61 


1.77 


1.51 


0.31 



tribution decreases and the centroid of the TLF* veloc- 
ity distribution increases as the reaction is increasingly 
damped. While for b>7 fm, a difference between the av- 
erage parallel velocity for IMFs and the PLF* (or TLF*) 
is discernible, for 4<b<7, the IMF distribution is cen- 
tered on (Vplf*)- The widths of the distributions are 
presented for completeness. No consistent trend of sig- 
nificance is evident in the extracted widths. 

We have also used the two gaussian fits previously 
described to extract the average multiplicity associated 
with the PLF*, TLF* and mid- velocity components at 
t=300 fm/c. The results are tabulated in Table|Il| With 
increasing centrality the multiplicities for each compo- 
nent increases although for the most central collisions, 



b< 7 fm the mid-velocity multiplicity seems to satu- 
rate at a value of sal. 5. From the peripheral collisions, 
8.5<b<10, to the mid-central collisions, 4 < b < 7, the 
average multiplicity of the PLF* and TLF* components 
increases by a factor «10 with an increase by ~25 for 
the mid- velocity component. The relative multiplicity 
of mid-velocity emission as compared to the total PLF* 
and TLF* emission increases from 0.26 for peripheral col- 
lisions to 0.44 for more central collisions. 

The velocity distributions of the IMFs are significantly 
altered by secondary decay. This influence is most evi- 
dent for the peripheral collisions where the shape of the 
primary distribution is nearly completely destroyed. Nat- 
urally, the magnitude of this secondary decay is particu- 
larly sensitive to the excitation predicted for the primary 
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FIG. 15: Left column: Parallel velocity distributions for the 
PLF* and TLF* (shaded), as well as IMFs (solid histogram) 
as a function of impact parameter. The dashed histograms 
correspond to a two gaussian fit as described in the text. 
Right column: Parallel velocity distributions of the PLF, 
TLF, and IMFs following Coulomb propagation and decay. 
The IMFs distributions have been normalized to the number 
of events. The PLF*, TLF*, PLF and TLF distributions have 
been scaled relative to the IMFs distributions. The arrows in- 
dicate the projectile and target velocities. 



fragments. The observed influence of secondary decay on 
the IMF velocity distribution indicates that the IMFs are 
significantly excited. 

For peripheral collisions, it has been experimentally 
observed that the emission pattern of a particles emitted 
by the PLF* manifests an anisotropic distribution |4{j[ . 
This observed anisotropy has been interpreted as the de- 
cay of a PLF* (and TLF*) initially deformed by the col- 
lision process. To investigate the extent to which such a 
physical picture is compatible with the AMD model, we 
have examined the invariant cross-section maps of a par- 
ticles. Depicted in Fig.^|is the dependence of the invari- 
ant cross-section map for a particles on impact parameter 
both at t=300 fm/c and following Coulomb propagation 
to infinite separation and sequential decay. At t=300 
fm/c (left column), for the most peripheral collisions the 
yield is peaked near the average PLF* and TLF* veloc- 
ities (indicated by arrows), though slightly toward the 
center-of-mass velocity. For these most peripheral col- 
lisions one observes that the primary a yield centered 
at mid-velocity is relatively small. With decreasing im- 



pact parameter, the primary alpha distributions associ- 
ated with the PLF* and TLF* move closer in velocity 
and increasingly overlap. 

Following secondary decay, the Coulomb circles evi- 
dent for a particles are even more striking than those for 
IMFs. This observation is consistent with the large multi- 
plicity of a particles that originate from the de-excitation 
of the PLF* and TLF* as compared to the early dynam- 
ical stage. The distinct emission pattern observed for 
8.5<b<10 fm is also observed for more central collisions 
although with increasing centrality the distinct nature 
of the semi-circles becomes less striking. The ridge of 
yield which is typically interpreted as Coulomb barrier 
emission becomes broader and its center moves increas- 
ingly toward Vn =0. These trends are consistent with 
the increased damping, excitation, and reduced size of 
the PLF* (TLF*) with increasing centrality. For even 
the most peripheral collisions, the pattern evident in the 
right column of Fig. 1161 is clearly not isotropic, favor- 
ing backward emission. As the sequential decay of the 
PLF* following t=300 fm/c is taken to be that of an iso- 
lated spherical nucleus without considering the influence 
of the external Coulomb field of the target on its decay 
[isj. it does not contribute to the predicted anisotropy. 
Within the model calculation, the observed anisotropy 
has two possible origins: Coulomb focusing of the a par- 
ticles present at t=300 fm/c and the a decay of IMFs 
which are emitted anisotropically. 

The anisotropic emission of a particles along the PLF* 
Coulomb ridge has recently been proposed to be related 
to the enhanced backward decay of the excited PLF* 
due to the nucleus-nucleus interaction |48j| . Displayed in 
Fig. El is the a particle yield along the Coulomb ridge 
for "peripheral" collisions, 5<b<10 fm, and "central" 
collisions, b<5 fm. Alpha particles were selected to be 
"Coulomb barrier" particles by restriction on their ve- 
locities, namely V a < 3.5 cm/ns in the PLF* frame. In 
both cases shown, the total a particle yield (solid his- 
togram) is not symmetric with respect to emission trans- 
verse to the PLF* direction, namely cos(6> Q )=0. Emission 
in the backward direction cos(0 a )<O is enhanced with re- 
spect to the forward direction. This enhancement is more 
pronounced for the central collisions. For the peripheral 
collisions, the emission yield for cos(9 a )= -1 is approxi- 
mately 1.7 times the yield emitted in the transverse di- 
rection. In contrast, the forward emission yield cos(# Q ) = 
+1 is approximately the same as the transverse yield. 
Comparison of the integrated yield with -l<cos(6> Q )<0, 
Ybackward{a), to 0<cos(6> a )<+l, Y forward (a) , reveals 
that backward emission is enhanced by ~19% as com- 
pared to forward emission. For more central collisions 
(bottom panel), comparison of the integrated yield re- 
veals that backward emission is enhanced by «39% as 
compared to forward emission. 

We have investigated the origin of this backward en- 
hancement, by examining the possible sources of a par- 
ticles. Alpha particles are "tagged" as being either a) 
"primary" , namely those originating at the time of clus- 
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FIG. 16: (Color online) Invariant cross-section for a particles in the COM frame. The arrows indicate the average parallel 
velocity of the PLF* and TLF*. The color scale indicates the yield on a logarithmic scale. In the left column the vertical scale 
of the bottom panel is scaled by a factor of six as compared to the upper panels. The right column is scaled by a factor of 
three with respect to the left column. 



ter recognition (t=300 fm/c), b) PLF* alphas or c) clus- 
ter alphas i.e. those that result from the secondary decay 
of primary IMFs. As expected, PLF* emission is essen- 
tially isotropic. It is evident in Fig. El that for both 
peripheral and central collisions, primary alphas (open 
triangles) on the PLF* Coulomb ridge are isotropic. Ev- 
idently the Coulomb focusing of primary a particles by 
the PLF* and TLF* does not contribute to the anisotropy 
observed in Fig. ^| The large backward enhancement 
observed for the total a particle yield is associated with 
the a particles that originate from the secondary decay of 
primary IMFs. Hence, it is the anisotropy of the primary 
IMFs that is responsible for the anisotropy of a particles 
associated with Coulomb barrier energies. Quantitative 
comparison of the the various components reveals that for 
peripheral collisions the ratio Y backward (a) /Y f orward (a) 
associated with the decay of primary clusters is »2.2, 
while the same ratio for PLF* emission or primary a 
emission is «0.9. In the case of more central collisions, 
the ratio Y backward (a) /Y f orward (a) associated with the 
decay of primary clusters is asl.7. The yield ratio for 
PLF* emission and primary a emission is as 1.0 in this 
impact parameter interval. 

The anisotropy observed for the decay of primary clus- 
ters is consistent with the emission pattern of IMFs as 



shown in Fig. This feeding of a particles to the PLF* 
Coulomb ridge from IMF secondary decay is also consis- 
tent with the fact that IMFs are produced excited, even 
for the most peripheral collisions. In fact, the average 
excitation energy of the IMFs is relatively independent 
of the impact parameter. For the most peripheral col- 
lisions, (E*/A) of the IMFs is typically 2.5 to 3 MeV 
with the higher values associated with IMFs produced 
around the center-of-mass velocity. With increasing cen- 
trality, (E*/A) becomes independent of the IMF velocity 
and reaches a typical value of 3 MeV. Such an excitation 
energy is in agreement with the excitation energy experi- 
mentally deduced for IMFs produced in central collisions 
|49|. Investigation of the width of the IMF excitation 
energy distribution reveals that it is large and approxi- 
mately independent of the impact parameter. 

In addition to the anisotropies predicted by the model, 
in reality the anisotropic emission pattern of a parti- 
cles can have additional origins. Although the PLF* is 
clearly deformed for t<300 fm/c (as shown in Fig. QJ, 
the statistical decay of the PLF* (and TLF*) is assumed 
to be isotropic. However, if the collision dynamics pref- 
erentially "prepares" the system in a configuration that 
favors emission toward the center-of-mass, the observed 
emission pattern will certainly be anisotropic. An ex- 
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FIG. 17: Angular distribution for a particles on the PLF* 
Coulomb ridge (V a < 3.5 cm/ns). All angles are in the PLF* 
frame. 



ample of such a favored configuration would be a di- 
nuclear configuration of the PLF* decaying into an IMF 
and residue with the IMF preferentially oriented toward 
mid-rapidity. If the di-nuclear configuration prepared lies 
outside the saddle point for such a system, then the ex- 
citation energy of the di-nuclcar configuration does not 
influence the decay probability and the decay is clearly 
non-statistical. However, if the di-nuclear configuration 
lies inside the saddle point, excitation energy does in- 
fluence the decay probability and the emission can be 
considered statistical. In this case, explicit treatment of 
the deformation within a statistical framework is neces- 
sary |38l ] . The observed anisotropy under such conditions 
will depend sensitively on the emission time relative to 
the rotational period of the di-nuclcar system. Of course 
such a schematic description of the binary decay of the 
PLF* could be extended to ternary and quaternary de- 
cays. It should also be noted that such short timescale 
emission when the nuclei are in proximity of each other 
and can also be influenced by tidal effects [37| . 



VI. SUMMARY AND CONCLUSIONS 

Using the AMD model, we have examined the dynam- 
ical phase of a heavy-ion collision at intermediate en- 
ergy. We have investigated how observables such as the 
size (Z), velocity, and excitation of the reaction products 



evolve during the early stages of the collision. The de- 
excitation of the initial reaction products is calculated 
with a statistical decay code and the survival of these 
initial observables is examined. 

We have investigated how the characteristics of the 
two large remnants in the reaction evolve with impact 
parameter. Both the (Z) and (V) of the PLF* and TLF* 
decrease smoothly as centrality increases up to an im- 
pact parameter of ps4 fm. As the centrality increases 
from the most peripheral collisions, the PLF*'s veloc- 
ity is increasingly damped from the projectile velocity. 
Concurrent with this damping, the width of the velocity 
distribution increases. Although the average velocity is 
largely unchanged by secondary decay, the width of the 
velocity distribution is typically increased by 10-40 %. 
For smaller impact parameters, b<3-4 fm, the average 
atomic number and velocity of the two reaction part- 
ners are independent of impact parameter. Associated 
with these changes in the size and velocity and size of 
the PLF* and TLF*, one also observes that the (E*/A) 
of the PLF* and TLF* increases as the impact parame- 
ter decreases from an initial value of 0.7-1.1 MeV upto 4 
MeV. The maximum excitation energy is attained for an 
impact parameter of ~6 fm. Smaller impact parameters 
do not result in larger values of (E*/A). These obser- 
vations suggest that the peripheral collisions on one side 
and the most central collisions on the other side corre- 
spond to different dynamics regime although simulated 
with the same ingredients. 

Peripheral collisions, as may be expected, exhibit a bi- 
nary nature with a strong memory of the entrance chan- 
nel. In such collisions, a transiently deformed PLF* and 
TLF* are recognizable as early as «100 fm/c after the 
collision. The deformation of these reaction products 
persists for a considerable time, t>300 fm/c. In addition 
to the PLF* and TLF*, nucleons, light charged parti- 
cles, and IMFs are also produced in the dynamical phase. 
The latter clusters, are preferentially located between the 
PLF* and TLF*. The Z distribution of particles with 
V|i>0 strongly favors asymmetric splits. The popula- 
tion of symmetric splits increases for mid-peripheral col- 
lisions reflecting an increase in the excitation energy of 
the PLF*. The excitation energy of the PLF* is strongly 
correlated with its velocity damping. Both the (E*/A) 
and (Vplf*) in this impact parameter range are found 
to be independent of the cluster recognition time. The 
correlation between the excitation energy of the PLF* 
and velocity damping is the same for the different clus- 
ter recognition times studied. The general inscnsitivity 
of (Vplf*) to cluster recognition time makes it a robust 
signal of the impact parameter. The excitation energy of 
the PLF* is slightly correlated with the excitation energy 
of the TLF* for early cluster recognition times. Even the 
small particle emission that occurs on short timescale is 
sufficient to destroy this correlation by t=300 fm/c. 

In contrast to peripheral collisions which exhibit a 
strong binary character, the most central collisions do 
not manifest as much memory of the entrance channel. 
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Such collisions are no longer dominated by two large 
fragments, namely the PLF* and TLF*. However, if we 
designate the largest fragment forward and backward of 
the center-of-mass as the PLF* and TLF*, their char- 
acteristics, (Z), (V) and (E*/A), are largely unchanged 
as b decreases for b<4 fm. Therefore, for the inner- 
most «10 % of the total cross-section, the maximum de- 
gree of excitation for such collisions is attained. This 
broad range of impact parameters associated with high 
excitation underscores the importance of considering the 
breakup of non-spherical geometries 50] . Moreover, for 
these small impact parameters, the quantitative charac- 
teristics of the PLF* and TLF*, (Z) and (E*/A), depend 
on the cluster recognition time. For early cluster recog- 
nition time (t=150 fm/c) an average excitation energy of 
6 MeV is reached, while a longer cluster recognition time 
(300 fm/c) results in (E*/A) of 4 MeV. This decrease 
in (E*/A) indicates a rapid de-excitation during the dy- 
namical stage, suggesting significant nucleon and cluster 
emission on a short timescale. 

Direct examination of the multiplicities of emitted 
particles reveals that the IMF multiplicity increases 
smoothly with increasing centrality and saturates for an 
impact parameter of «3 fm. At t=300 fm/c, the average 
IMF multiplicity reaches unity for an impact parameter 
of «9 fm. At all impact parameters and for t> 150 fm/c, 
the IMF emission rate decreases monotonically with in- 
creasing cluster recognition time. The velocity distribu- 
tion of the produced IMFs at t=300 fm/c is bimodal and 
reveals preferential emission from the PLF* and TLF* to- 
wards the center-of-mass. The emission pattern of a par- 
ticles at t=oo also exhibits a distinct preferential emis- 
sion towards the center-of-mass. This anisotropy how- 
ever, is not due to the anisotropic emission of primary a 
particles or evaporation from the PLF* and TLF*, but 
arises from the secondary decay of anisotropically emit- 
ted primary IMFs. The multiplicity of neutrons saturates 



for mid-central collisions, making neutron multiplicity a 
poor selector of central collisions. The geometric cross- 
section combined together with the saturation of excita- 
tion energy for mid-peripheral collisions may explain the 
observed persistence of binary collisions at intermediate 
energies even when the largest neutron multiplicities are 
selected |47j . 

The large excitation energy reached in the colli- 
sion leads to rapid particle emission on the dynamical 
timescale. However, the present treatment of the short 
timescale decay involves several simplifications. The role 
of deformation in the decay is neglected as are both 
Coulomb and nuclear proximity effects. In addition the 
excitation energy is calculated relative to clusters that 
are at the ground-state both in shape and density. These 
simplifications may have a non-negligible impact on the 
characteristics of the fragmenting system. This rapid de- 
excitation emphasizes the need for a hybrid statistical- 
dynamical model that considers in a more realistic man- 
ner the statistical decay of the transiently deformed nu- 
clei from times as short as 100 fm/c. Development of 
such a hybrid model would represent a new and poten- 
tially powerful tool in understanding the dynamics of in- 
termediate energy heavy-ion collisions, as well as cluster 
formation on short timescales. 
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